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I discuss recent progress in understanding the high-energy evolution in QCD, which points 
towards a remarkable correspondence with the reaction-diffusion probfem of statistical physics. 



Recently, there has been significant progress in our understanding of QCD at high energy, 
based on the observations that (i) the gluon number fluctuations play an important role in 
the evolution towards saturation and the unitarity limit ^'^ and (ii) the QCD evolution in the 
presence of fluctuations and saturation is in the same universality class as a series of problems 
\ in statistical physics, the prototype of which being the 'reaction-diffusion' problem These 

' observations have developed into a rich correspondence between high-energy QCD and modern 

. problems in statistical physics, which relates topics of current research in both fields, and which 

\ has already allowed us to deduce some insightful results in QCD by properly translating the 

' corresponding results from statistical physics 

To put such theoretical developments into a specific physical context, let us consider 7*- 



proton deep inelastic scattering (DIS) at high energy, or small Bjorken-x. We shall view this 
process in a special frame in which most of the total energy is carried by the proton, whose 
wavefunction is therefore highly evolved, while the virtual photon has just enough energy to 
, dissociate long before splitting into a quark-antiquark pair in a colorless state (a 'color dipole'), 

which then scatters off the gluon distribution in the proton. The transverse size r of the dipole is 
controlled by the virtuality of 7* (roughly, ~ so for ^> Aqqj-, one can treat the 

dipole scattering in perturbation theory. But for sufficiently small x, even such a small dipole 
5^ \ can see a high-density gluonic system, and thus undergo strong scattering. 

Specifically, the small-x gluons to which couple the projectile form a color glass condensate^ 
(CGC), i.e., a multigluonic state which is characterized by high quantum occupancy, of order 
l/os, for transverse momenta k± below the saturation momentum Qs{x), but which becomes 
rapidly dilute when increasing k± above Qs- The saturation scale rises very fast with the energy, 
Q'^{x) x"'^, and is the fundamental scale in QCD at high energy. In particular, the external 
dipole is strongly absorbed provided its size r is large on the scale set by 1/Qs, whereas for 
r <C i/Qs one rather has weak scattering, or 'color transparency'. 

In turn, the small-x gluons are produced through quantum evolution, i.e., through radiation 
from color sources (typically, other gluons) with larger values of x, whose internal dynamics is 
'frozen' by Lorentz time dilation. Let r = In 1/x denote the rapidity ; it takes, roughly, a rapidity 
interval At ~ 1/ces to emit one small-x gluon; thus, in the high energy regime where a^r ^ 1, 
the dipole meets with well developed gluon cascades, as illustrated in Fig. 1. Three types of 
processes can be distinguished in Fig. 1, which for more clarity are singled out in Fig. 2. 

The first process. Fig. 2. a, represents one step in the standard BFKL evolution^; by iter- 
ating this step, one generates gluon ladders which are resummed in the solution to the BFKL 
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Figure 1: An instantaneous gluon conSguration in the proton wavefunction as 'seen' in DIS at small x. 




equation^. However, by itself, the latter is well known to suffer from conceptual difficulties in 
the high energy limit : (i) The BFKL estimate for the dipole scattering amplitude ^^(r) grows 
exponentially with r (i.e., like a power of the energy), and thus eventually violates the unitarity 
bound < 1. (The upper limit T,- = 1 corresponds to the black disk limit, in which the 

dipole is totally absorbed by the target.) (ii) The BFKL ladder is not protected from devia- 
tions towards the non-perturbative domain at low transverse momenta ^ ^qcd ('infrared 
diffusion'). With increasing energy, the BFKL solution receives larger and larger contributions 
from such soft intermediate gluons, and thus becomes unreliable. 

These 'small-x problems' of the BFKL equation are both cured by the gluon recombination 
processes {n — > 2) illustrated in Fig. 2.b which are important at high energy, when the gluon 
density in the target is large, and lead to gluon saturation and the formation of the CGC. Such 
processes are included in the Balitsky-JIMWLK equation^, a non-linear, functional, general- 
ization of the BFKL evolution which describes the approach towards gluon saturation in the 
target and preserves the unitarity bound in the evolution of the scattering amplitudes. 

However, the Balitsky-JIMWLK equation misses the process in Fig. 2.c — the 2 — > n 
gluon splitting — which describes the bremsstrahlung of additional small-x gluons in one step 
of the evolution. By itself, this process is important in the dilute regime, where it leads to the 
construction of higher-point gluon correlation functions from the dominant 2-point function. 
But once generated, the ra-point functions with n > 2 are rapidly amplified by the subsequent 
BFKL evolution (the faster the larger is n) and then they play an important role in the non- 
linear dynamics leading to saturation. Thus, such splitting processes are in fact important for 
the evolution towards high gluon density, as originally observed in numerical simulations ^ of 
Mueller's 'color dipole picture'^, and more recently reiterated in the analysis in Refs. ^'2''^. 

Equations including both merging and splitting in the limit where the number of colors Nc 
is large have recently became available^ (see also Refs. ^'^), but their general solutions have 
not yet been investigated (except under some additional approximations^' ^'^). Still, as we shall 



argue now, the asymptotic behaviour of the corresponding solutions — where by 'asymptotic' 
we mean both the high-energy limit r ^ oo and the weak coupling limit ^ — can be 
a priori deduced from universality considerations, by exploiting the correspondence between 
high-energy QCD and the reaction-diffusion problem of statistical physics^. 

To that aim, it is convenient to rely on the event-by-event description ^ of the scattering 
between the external dipolc and the hadronic target (cf. Fig. 1) and to use the \argc-Nc 
approximation to replace the gluons in the target wavefunction by color dipoles ^. Then, the 
dipole-target scattering amplitude corresponding to a given event can be estimated as 



where is the scattering amplitude between two dipoles with comparable sizes and nearby 
impact parameters, and frir, h) is the occupation number for target dipoles with size r at impact 
parameter 5, and is a discrete quantity: / = 0, 1,2, . . . . Thus, in a given event, the scattering 
amplitude is a multiple integer of a^. 

In this dipole language, the 2 — 4 gluon splitting depicted in Fig. 2.c is tantamount to 1 — > 2 
dipole splitting, and generates fluctuations in the dipole occupation number and hence in the 
scattering amplitude. Thus, the evolution of the amplitude TT-(r, h) with increasing r represents 
a stochastic process characterized by an expectation value (T(r, h))r — oil {f{r, b))r, and also by 
fluctuations ST ~ oi^Sf ~ \/ct^T (we have used the fact that 5f ~ \ff for fluctuations in the 
particle number). Clearly, such fluctuations are relatively important (in the sense that 5T ^ T) 
only in the very dilute regime where (/) ^ 1, or (T) ^ a^. 

Eq. (1) applies so long as the scattering is weak, T ^ 1, but by extrapolation it shows that the 
unitarity corrections are expected to be important when the dipole occupation factor becomes 
of order l/a^. Consider first the formal limit a1 0, in which the maximal occupation number 

^ ^/oil becomes arbitrarily large. Then one can neglect the particle number fluctuations 
and follow the evolution of the scattering amplitude in the mean field approximation (MFA). 
This is described by the Balitsky-Kovchegov equation ^, a non-linear version of the BFKL 
equation which, as shown in Rcf. ^, lies in the same universality class as the Fishcr-Kolmogorov- 
Petrovsky-Piscounov (FKPP) equation (the MFA for the reaction-diflFusion process and related 
phenomena in biology, chemistry, astrophysics, etc; see^^ for recent reviews and more references). 
The FKPP equation reads, schematically. 



in notations appropriate for the dipole scattering problem: T(/9, r) = {T{r))r and p = ln(rQ/r^), 
with ro a scale introduced by the initial conditions at low energy. Note that weak scattering 
(T ^ 1) corresponds to small dipole sizes (r <C 1/Qs), and thus to large values of p. The 
three terms on the r.h.s. of Eq. (2) describe, respectively, diffusion, growth and recombination. 
The flrst two among them represent (an approximate version of) the BFKL dynamics, while 
the latter is the non-linear term which describes multiple scattering and thus ensures that the 
evolution is consistent with the unitarity bound T < 1. 

Speciflcally, the solution Tt{p) to Eq. (2) is a front which interpolates between two fixed 
points : the saturation fixed point T = latp^—oo and the unstable fixed point T = at 
p ^ oo (see Fig. 3). The position of the front, which marks the transition between strong 
scattering (T ~ 1) and, respectively, weak scattering (T <C 1), defines the saturation scale: 
Ps{t) = \n{rQQl{T)) . With increasing r, the front moves towards larger values of /?, as illustrated 
in Fig. 3. Note that the dominant mechanism for propagation is the BFKL growth in the tail of 
the distribution at large p : the front is pulled by the rapid growth of a small perturbation around 
the unstable state. In view of that, the velocity of the front A = dpg/dr is fully determined by 
the linearized version of Eq. (2), which describes the dynamics in the tail. Specifically, by solving 



r,(r,6) ~ alU{r,b), 



(1) 



drT{p,T) = dlT{p,T) + T(p,r) - T\p,t) 



(2) 
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Figure 3: Evolution of the continuum front of the Balitsky-Kovchegov equation with increasing r. 



the BFKL equation one finds ^'^^'^^ that, for p > Ps{t) and sufficiently large r, 

Trip) - e'^'^'^e-^P = e-^^^-'^^M) , Ps{t) = ca^r, (3) 

where = agNc/Tr, 7 = 0.63.., and c = uj/j = 4.88... From Eq. (3) one can immediately 
identify the velocity of the front in the MFA as Aq = cctg. Since Q'^{t) ~ Qq e^""^, it is furthermore 
clear that Aq plays also the role of the saturation exponent (here, in the MFA). 

What is the validity of the mean field approximation ? We have earlier argued that the 
gluon splitting processes (cf. Fig. 2.c) responsible for dipole number fluctuations should play 
an important role in the dilute regime. This is further supported by the above considerations 
on the pulled nature of the front: Since the propagation of the front is driven by the dynamics 
in its tail where the fluctuations are a priori important, the front properties should be strongly 
sensitive to fluctuations. This is indeed known to be the case for the corresponding problem in 
statistical physics ^^'j as it can be understood from the following, qualitative, argument: 

Consider a particular realization of the stochastic evolution of the target, and the corre- 
sponding scattering amplitude which is discrete (in steps of a^). Because of discreteness, the 
microscopic front looks like a histogram and thus is necessarily compact : for any r, there is 
only a finite number of bins in p ahead of ps where Tj- is non-zero (see Fig. 4). This property 
has important consequences for the propagation of the front. In the empty bins on the right of 
the tip of the front, the local, BFKL-like, growth is not possible anymore (as this would require 
a seed). Thus, the only way for the front to progress there is via diffusion, i.e., via radiation 
from the occupied bins at p < pup (compare in that respect Figs. 3 and 4). But since diffusion 
is less effective than the local growth, we expect the velocity of the microscopic front (i.e., the 
saturation exponent) to be reduced as compared to the respective prediction of the MFA. 

To obtain an estimate for this effect, we shall rely again on the universality of the asymp- 
totic (r — > 00 and = l/a^ ^ 1) behaviour^. Namely, from the experience with the reaction- 
diffusion process and related problems in statistical physics one knows indeed that the 
dominant behaviour for large evolution 'time' and large (but finite) occupancy A ^ 1 is inde- 
pendent of the details of the microscopic dynamics, and thus is the same for all the processes 
whose mean field limit {N — > 00) is governed by the FKPP equation (2). In particular, the 
dominant contribution to the correction Aat — Aq to the front velocity is known to be universal, 
and can be obtained through the following, intuitive, argument, due to Brunet and Derrida^^: 

For a given microscopic front and A ^ 1, the MFA should work reasonably well everywhere 
except in the vicinity of the tip of the front, where the occupation number / becomes of order 
one (corresponding to T ~ in the QCD problem) and thus the linear growth term becomes 
ineffective. Accordingly, Brunet and Derrida suggested a modified version of the FKPP equation 
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Figure 4: Evolution of the discrete front of a microscopic event with increasing rapidity r. The small blobs are 
meant to represent the elementary quanta of T in a microscopic event. 



(2) in which the 'BFKL-hke' growth term is switched off when T < a 

drT{p,T) = d^T + G(r-a2)r(i-r). 
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(4) 



By solving this equation in the hnear regime, they have obtained the first correction to the front 
velocity as compared to the MFA (in notations adapted to QCD; see Ref. ^ for details): 
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(5) 



where the numbers c « 4.88 and n ~ 150 are fully determined by the linear (BFKL) equation. 
In QCD, the same result has been first obtained through a different but related argument by 
Mueller and Shoshi ^ . Note the extremely slow convergence of this result towards its mean field 
limit: the corrective term vanishes only logarithmically with decreasing = a^, rather than 
the power-like suppression usually found for the effects of fluctuations. This is related to the 
high sensitivity of the pulled fronts to fluctuations, as alluded to above. Such a slow convergence, 
together with the relatively large value of the numerical factor k, make that the above estimate 
for A, although exact for asymptotically small a^, is in fact useless for practical applications. To 
understand the subasymptotic corrections and, more generally, the behaviour of the saturation 
momentum and of the scattering amplitudes for realistic values of r and , one needs to solve 
the actual evolution equations of QCD^'^'^, a program which is currently under way^'^. 
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